Spectral Analysis and Identification of Noises in Quantum Systems 
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In quantum information processing, knowledge of the noise in the system is crucial for high- 
precision manipulation and tomography of coherent quantum operations. Existing strategies for 
identifying this noise require the use of additional quantum devices or control pulses. We present a 
noise-identification method directly based on the system's non-Markovian response of an ensemble 
measurement to the noise. The noise spectrum is identified by reversing the response relationship 
in the frequency domain. For illustration, the method is applied to superconducting charge qubits, 
but it is equally applicable to any type of qubits. We find that the identification strategy recovers 
the well-known Fermi's golden rule under the lowest-order perturbation approximation, which corre- 
sponds to the Markovian limit when the measurement time is much longer than the noise correlation 
time. Beyond such approximation, it is possible to further improve the precision at the so-called 
optimal point by incorporating the transient response data in the non-Markovian regime. This 
method is verified with experimental data from coherent oscillations in a superconducting charge 
qubit. 

PACS numbers: 03.67.-a, 42.50.Lc, 85.25.Dq 



I. INTRODUCTION 



Artificial atoms fabricated with solid-state devices 
(e.g., superconducting qubits, quantum dots, and NV 
centers) are promising for future quantum information 
processors (e.g., pHEJ an d references therein). To con- 
struct more stable and robust quantum circuits, much ef- 
fort is required to overcome complex decoherence brought 
by their "dirty" environment. Higher performance can 
only be achieved by taking control actions (e.g., the dy- 
namical decoupling or environment modulation (Tol — 
[lH) to fight against decoherence. To characterize and 
control quantum systems, it is crucial to obtain a dy- 
namical model of the quantum system according to the 
measurement result of the system. For closed systems, 
this calls for the identification of the Hamiltonian [l4| - [l6| . 
Otherwise, when the system is open, the dissipation ef- 
fect also needs to be specified corresponding to the noises 
coupled to the system. Hence, the information acquisi- 
tion of the noises is critical for modeling open quantum 
system dynamics. The knowledge of noise is also useful 
for improving the fidelity of state and process tomog- 
raphy from a series of designed quantum measurements 
jl7H2l| . Since the number of required measurements for 
process tomography increases quickly with the system 
dimensionality, a good model for the dynamics of the 
measured system will help reduce the measurement cost 
based on the prediction of the actual process matrix. 

Physically, a stationary noise in a quantum system 
is mainly characterized by its correlation spectrum [22| . 
In many solid-state systems, the spectrum has nontriv- 
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FIG. 1: (color online) Schematics for a general identification 
procedure. The quantum dynamics is properly excited so that 
the noise n can have a simple relation with the measurement 
data m, which is then reversed to numerically identify the 
noise properties. 



ial structures (e.g., Lorentzian shape, Ohmnic, 1/f noise 
[23j). which corresponds to colored noise. Figure Q] shows 
the general procedure for identifying an unknown noise 
"n" acting on a quantum system. The expectation value 
m of some observable is measured. If one knows how the 
noise n affects m via a functional mapping m = S^ip) 
[e.g., Eq. (|18p in Sec lIIIj . then this mapping can be re- 
versed (if reversible) to identify the noise from the mea- 
surement data, i.e., n = I 5^~ 1 (m) [e.g., Eq. (fTTj)) in 
Sec lIIIj . Hence, the identification precision largely de- 
pends on how well the model m = S^(vl) describes the 
real noise-drive quantum dynamics. 

Generally, the noise causes a random frequency shift 
and dissipation effects in the system dynamics, which can 
be described by a generalized master equation 0, 
or a spectral overlap function with the noise correla- 
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tion spectrum [12J, |24[ . Such models were used to derive 
l2r| the widely used identification strategy based on 
Fermi's golden rule, which claims that the state transition 
rate of a two-level system is proportional to the intensity 
of the noise spectrum at the resonating frequency. This 
makes it possible to experimentally anal yze the noise by 
using a frequency-tunable qubit [3j, |25U28| . Moreover, 
in Refs.pS. |3G |. the dephasing noise spectrum is identi- 
fied from the direct measurement of the asymptotic de- 
cay of the nuclear spin qubit in NMR experiments. The 
strategy requires designed 7r-pulse sequences to cancel the 
dominant influence of the lower-frequency component of 
the noise spectrum (3j [HJ 

The above identification schemes essentially exploit 
only the system's long-time (compared to the noise cor- 
relation time) response to the noise [25j . while the non- 
Mar kovian transient dynamics is not considered. In this 
paper, we introduce a new identification method that 
incorporates the non-Markovian noise-driven dynamics. 
The model adopted for identification is based on the Born 
approximation, which has been used to compute non- 
Markovian dynamics driven by known colored noises (33l — 
[36j . Under the Laplace transform, the model appears in 
a purely algebraic form that will simplify the analysis for 
the purpose of noise identification. 

This paper is organized as follows. In Section [TTl a 
frequency-domain model is derived for analyzing the re- 
sponse of any ensemble measurement to the noise in a 
quantum system. In Section IIII[ the model is applied 
to a qubit system and is shown to be consistent with 
Fermi's golden rule. Then, from the model we propose 
a new identification strategy at the optimal point and 
demonstrate it with experimental data. Finally, Section 
IIVI presents our conclusions. 



II. FREQUENCY-DOMAIN MODELING OF 
OPEN QUANTUM SYSTEMS 

In this section, we will introduce the generalized mas- 
ter equation based on the Born approximation, and 
Laplace transform the equation to the frequency-domain 
for the following identification analysis. Although such 
model has been obtained in different ways in the litera- 
ture (34l - [3l| , we will express it in a new compact form to 
highlight the influences of the noise on the qubit dynam- 
ics. This will facilitate the following study on identifying 
the noise spectrum from the measurement data of a phys- 
ical observable. 

Consider an iV- level quantum system that is coupled 
to some unknown environment. Suppose that the total 
Hamiltonian is 



H Q + H 1 ^C{t), 



where Hq is the internal Hamiltonian of the system. The 
environment operator C(t) coupled to the system's ob- 
servable H\ is responsible for the noise. For convenience, 
it has been rotated to the interaction picture. Assume 



that the environment is initially at an equilibrium state 
p B and the coupling strength between the system and 
the environment is so weak that the joint evolution with 
the system hardly alters the state of the environment. 
This qualifies the Born approximation [34j . under which 
the net system evolution (after averaging out the envi- 
ronment) can be described by the follo wing generalized 
master equation (with H=l) pH[l2l [3H37j: 



p{t) = ~i [H 0>P (t)]- ^dt'{$(t,t') 
H 1 p{t')e i( - t - t ^ H °'] +H.c.}, 



-i(t-t')H . 



(1) 



where p(t) is the density matrix of the system. 

The integral term in the above equation captures the 
noise-driven quantum dynamics, where the environmen- 
tal noise enters through the temporal noise correlation 
function 

t') = Tr [C(t)C(t')p° B ] = T(t, t') + iSl(t, t'), (2) 

where T(t, t 1 ) and f2(f , £') are the real and imaginary parts 
of $(t,i'), respectively. Because the environment is ini- 
tially prepared at an equilibrium state p B , the correlation 
function only depends on the time difference, t—t'. Thus, 
in the following, $(•) [as well as T(-) and will t> e al- 
ways taken as a single- variable function of r = t — t' . 
The system behaves in a non-Markovian manner when 
the time is comparable with the noise correlation time. 

Here, V(t, t') and tt(t,t') in the correlation func- 
tion $(£,£') correspond to the commutative and non- 
commutative parts of C(t), respectively, as follows 

r(r) = ±Tr{[C(T)C(0)+C(0)C(T)]p° B }, 

n(r) = ^Tv{[C(r)C(0)-C(0)C(r)}p%}. 

In particular, when the environment is at zero temper- 
ature, T(t) corresponds to the average dissipation rate, 
while fl(t) corresponds to the Lamb shift. 

To facilitate the following analysis, we first transform 
Eq. ([l} for the density matrix p(t) into a vector form. 
Choose {M = N^In, Mi,- ■■ ,M N *_i} as an orthonor- 
mal basis of the space of N x N Hermitian matrices, 
where each Mi is Hermitian and Tr (MjMj) = 6ij for any 



< i < j < N 2 
be spanned as 



1. Then, the density matrix p(t) can 



p(t) = v (t)M 



«jv2_i(i)M 



N 2 -l, 



where Vi(t) — Tr[p(t)M;]. The real vector v(t) = 
[vo(t), • • ■ , v N 2_ 1 (t)] T is called the augmented Bloch rep- 
resentation of p{t) (see Section lTTTl for the example of two- 
level systems). In this way, Eq. |T]) can be translated into 
the following form: 



v(t) = C v(t) + Ci j dt' |r(t - i')e (t " t,)Co /;i 

+n(i-t')e (t ~'' )£o Ct}v(t'), (3) 
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where Ck and , k = 0,1, are the matrix represen- 
tations of the commutator operation [— and the 
anti-commutator operation {-fffc, •}, respectively. 

Eq. ((3]) can be naturally Laplace transformed to the 
frequency domain as only a linear term and a convolution 
term of v(t) are involved on the right hand side. This 
gives 

sv(s) -v° = C v(s) + C x K{s)v(s), (4) 

where v° is the Bloch vector corresponding to the initial 
density matrix. The matrix s-function 

K(a) = ^[r(r)e r£o ] Ci +^p{T)e TCo ] C+ 

= r(sl - £ )A + Q(sl - C )Cf (5) 

characterizes the action of the noise on the system evo- 
lution, where Jz?[-] denotes the Laplace transform and I 
is the identity operator on the space of Bloch vectors. 
Here the property Jzf [c(i)e s °*] = c(s — sq) of the Laplace 
transform of a scalar function c(t) has been extended to 
matrix- valued functions of L(sl — Co) and 0(sl — Co) (see 
Appendix A for details). 

Let us now denote by the matrix s-function TZ(s) = 
(si — £o) _1 the resolvent operator for the unperturbed 
system evolution. One can immediately derive from 
Eq. (0| the following frequency domain model: 

v( s )=[lZ- 1 ( S )-£ 1 JC( S )}- 1 v 1 (6) 

as well as that of the time derivative j(t) = v(t) (for 
calculating the transition rates) of the state vector: 

7(s) = sv(s) — v 

= [C + d)C(s)) ■ [n-\s) - A/C(s)] _1 v°. (7) 

The noise affects the system evolution through the opera- 
tor /C(s) as a perturbation to the coherent evolution IZ(s). 
Moreover, the fact that /C(s) is a function of si— Cq shows 
that the noise-driven dynamics is mainly determined by 
the properties of the noise near the system's frequencies. 

With application to actual experiments, Eqs. ([5]) and 
([7]) link the noise spectrum and the measurement data of 
some observable of the system, whose expectation value 
is always a linear function of v(t) or "f(t). The expression 
of such functional dependence is generally very compli- 
cated. However, as will be seen in the next section, it 
can be simplified under proper selection of the parame- 
ters and the system's initial state, so that useful identi- 
fication formulas can be derived. 



III. NOISE IDENTIFICATION IN A 
SUPERCONDUCTING CHARGE QUBIT 

This section will apply the above model to the identifi- 
cation of the noise spectrum in an open quantum system. 
As a physical example, we study two-level systems imple- 
mented by a superconducting charge qubit [3814421 ] . The 



method developed here can be easily extended to other 
physical implementations. 

As shown in Fig. [21 the superconducting charge qubit 
is encoded by the charge stat es |0) a nd |1) differing by 
one Cooper pair in the box [H, mSHl ■ The box is biased 
by the control gate voltage U via the capacitor C. Let Ej 
be the tunneling energy of the Joscphson junction, Eq be 
the Column energy of the island; and E c \ = Ec(l — 2n ff ) 
be the electrostatic energy that is linear in the dimen- 
sionless gate voltage n g (proportional to U in Fig. 
Then under the charge state basis {|0), |1)}, the qubit's 
internal Hamiltonian reads 

fftot = Eja x + E cl a z + a z ® C(t), 

where the bath operator C(t) is responsible for the charge 
noise (assumed to be dominant in experiments). It is 
commonly recognized that the noise is mainly induced 
by fluctuating background charges in the substrate or 
on the surface |46l448j . but in the model it can consist 
of any sources as we are only concerned with the over- 
all correlation properties. The average charge number, 
which corresponds to the expectation value of a z opera- 
tor, is read out by electrostatically coupling the qubit to 
a radio- frequency single-electron transistor (SET). 
The density matrix under the Pauli matrix basis 

p{t) = voit^+v^+v^+vsit)^ 

is vectorized as a four-dimensional augmented Bloch vec- 
tor v(t) = [v (t),v 1 ( t),v 2 (t),v 3 (t)} T with v (t) = 1. 
Moreover, let A = yj E 2 } + E^ } be the energy gap and 
9 = arctan(£'j / E c \) be the bias angle, we can rewrite the 
total Hamiltonian as 

fftot = A ere + <?z ® C{t), (8) 

where a 9 = a z cos 6 + a x sin 0. Next, we assume that 
the noise correlation time is much shorter than the de- 
cay time of the excited state, which validates the Born 
approximation [lj|, [3?J ■ The spectral relationship will be 
studied between the measurement result and the noise 
spectrum show under different bias angles and initial 
states, from which noise identification strategies can be 
designed. 

A. Fermi's Golden Rule 

The Fermi's golden rule has been broadly applied in 
spectroscopy techniques. In the application to supercon- 
ducting charge qubit systems [401 ) . the qubit is initially 
prepared at the ground state (or the excited state) under 
some bias voltage, which can be expressed as: 

6 6 6 6 

1+) =cos-|0)+sin-|l), I-) =-sin-|0)+cos-|l). 

The probability for the qubit to stay in the ground state 
(or the excited state) is measured at sampled time in- 
stances, from which the excitation (or decay) rate 7f(t) 
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FIG. 2: Schematic diagram for the circuit of a supercon- 
ducting charge qubit. The qubit is encoded by the num- 
ber of Cooper pairs in the box inside the dashed rectangle 
(with Josephson energy Ej) that is biased by the gate volt- 
age U (proportional to a dimensionless number n g of charges) 
through the capacitor (with capacitance C). 



[or 7j_ (t)J is measured by interpolating the time- variant 
curves. 

To model this process, it is convenient to express the 
total Hamiltonian (|8]) in the eigenbasis {|+), | — )}, which 
becomes: 



i? tot = Acr z + ao <g> C(t), 



(9) 



Both 7f(s) or 74. (s) can be calculated from the z- 
componcnt 



73(a) = - [sv 3 (s) 



of 7(s), where v® = 1 for 7f(s) and 



,0 - 



-1 for 74.(3), 



respectively. Their expressions under any bias angle are 
given by (|B2[) in Appendix [B] Because the noise oper- 
ator C(t) is weak under the Born approximation, it is 
reasonable to approximate Eq. (|B2|) by keeping only the 
lowest-order terms of T(s) and fi(s). This gives 



7t( s ) 
74 ( s ) 



sin 



[r + ( s ) + o_( s )], (10) 

[r + (a)-0_(a)], (11) 



where r+(s) = JSf [r(i) cos A<] and fi_(s) = 
Jzf [O(t) sin At] are the Laplace transforms of the mod- 
ulated noise correlation functions with a carrier wave 
whose frequency is A. Such approximation holds when 
the noise correlation time is much shorter than the decay 
time of the excited state. One can immediately see that 
the transition rates from |+) and |— ) are differentiated 
by the noise term f2_(s). When the time is sufficiently 
long (far greater than the noise correlation), the transi- 
tion rate from — ) to |+) approaches a stationary value, 
which can be calculated by taking the limit t — > 00 in 
Eqs. (|TU|) and (fTTj) . In the the frequency domain, these 



become 

7t = !ims7t( s ) 

s— >-0 

w sin 2 (9[r + (0) + fi_(0)] = $ FT (-A)sin 2 6>, (12) 
71 = lims7l(s) 

s— >0 

« sin 2 9 [r+ (0) - (0)] =$ft( A) sin 2 9, (13) 

where $ft(A) is the Fourier transform of the correla- 
tion function $ (£) 1491 . This exactly recovers the formu- 
las adopted in [26|, |4(| according to the Fermi's golden 
rule, which holds when the measurement time t is much 
smaller than the inverse of the stationary decay rate, and 
far greater than the autocorrelation time of the noise. 

The Fermi's golden rule shows that the spectral den- 
sity of a noise at some frequency A is proportional to the 
transition rate of a resonating two-level system. Thus, 
a transition-frequency tunable qubit can be used as a 
spectrum analyzer. In the charge qubit system, the tran- 
sition frequency A = Ej sin -1 9 is tuned by the bias an- 
gle 9. The identification scheme based on the Fermi's 
golden rule is operationally convenient, but the preci- 
sion may be limited due to the use of perturbation and 
asymptotic approximations, which equivalcntly requires 
that the system behaves approximately Markovian at the 
time of measurement (i.e., the environment has no mem- 
ory effects on the system's evolution). 



B. Noise identification at the optimal point 

The identification formula based on the Fermi's golden 
rule is obtained by measuring the decay rate under var- 
ious bias angles 9. Actually, it is possible to find better 
identification schemes from the full measurement data 
under a fixed bias angle. For example, the identification 
formula Eq. (|B2[) can be greatly simplified at 9 = ir/2 
(i.e., the so-called optimal point in superconducting qubit 
systems) as below 

7t( s ) = /x , . 74 (s) = -pi / \ 1 ■ (14) 



T + (s) + s 



r+( S )+.s 



The transition rates differ only by H_(s) in the numera- 
tors, which has its physical origin from the non-uniform 
thermal distribution of populations on the qubit levels. 
Eq. (|14[) can be easily reversed to obtain an identification 
formula: 



r+(*) = 

0_(s) = 



s[7t(s) +74 00] 
2-[7t00 + 7400]' 

s[7tOQ -74W] 
2 -[7t 00-7400]' 



(15) 
(16) 



for r_|_(s) and f2_(s), respectively, corresponding to the 
modulated signals r(t)cosAt and r2(t)sinAt, respec- 
tively. Because no asymptotic approximation is taken 
here, the identification is expected to be more precise 
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and more economic because only measurement data at 
the optimal point are required. 

The above formula can be used to identify the modu- 
lated signals r(i)cosAi [or f2(i) sin Ai]. However, from 
them one cannot obtain the full correlation function be- 
cause the values of T(t) [or f2(i)] at the time instances 
t = (k + h) -J (or t = ^r), where k is an arbitrary 
nonnegative integer, cannot be determined. Also, the 
presence of high-frequency measurement noise may cause 
large error in the calculation of 7f(i) or 74. (t). Neverthe- 
less, this idea can be extended to the following different 
experimental setup j38| , with which at least the spectrum 
of T(t) can be identified at the optimal point. 

As shown in Fig. Ela), the qubit is prepared at the 
zero-charge state |0) and released for free evolution under 
some bias angle 9. The average charge number is then 
measured at sampled time instants [see Figj3|b)], from 
which coherent oscillations can be observed [38[. The 
presence of charge noise C(t) will damp the oscillation 
via. 

Because the average charge number is measured, it is 
more convenient to use the original total Hamiltonian © 
in the charge states basis {|0}, |1)}. Correspondingly, 
the average charge number Qe(t) is calculated from the 
z-component of the Bloch vector: 



i(t) = \[l-va(!)] -> Qe{s) = \ 



- - V3(s) 

s 



(17) 



In AppcndixICl the functional dependence of Qg (t) on the 
noise spectrum is derived [see Eq. (|C2[) ] . At the optimal 
point 6 = I (corresponding to E e \ = 0), the formula 
becomes 



?tt/2(s) 



A 2 



2s [s 2 + sT{s) + A 2 ] ; 



(18) 



where only T(s) is present but Q(s) disappears. This can 
be used to separate the identification of T(t) and £l(t). By 
reversing Eq. (|T5|) at s = iui, we obtain an identification 
formula for T(t) 



CL 

03 

bo 
>■ 
no 

1_ 

c 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 
1 



o 0.8 




1000 2000 
t(ps) 



3000 



FIG. 3: (color online) Schematic diagram for identifying the 
noise T(t). (a) The energy gap versus the bias angle, where 
both quantities are defined below Eq. ([9]). The charge qubit 
is initially prepared in the zero-charge state, corresponding 
to the ground state far from the optimal point. Then, a gate 
pulse is applied to induce coherent oscillations at some bias 
angle 6 (e.g., the optimal point = -|), which is determined 
by the height of the pulse. The pulse persists for a time t, 
after which the average charge number is measured, (b) The 
coherent oscillation curve (with data collected from NEC and 
RIKEN) of the average charge number Q- K /2{t) versus the time 
t, which is measured at the optimal point 9 — tt/2. 



T(iw) 



2uj 2 Q 7r/2 {iuj 



a |, (19) 



Qir/2(t) ~~ 0-5 of the measurement data Q 7r /2(t)i which 
leads to the following identification formula 



where luq = Ej/H; and is obtained by Laplace 

transforming the measured coherent oscillation curve in 
the time domain. From (|19|) . the Fourier transform of 
the noise spectrum, which is more often used in practice, 
can be obtained as 



r FT (w) = 2 Rer(iw) 



-cjqRc 



1 



u 2 Q^ /2 {iuj) 



(20) 



where the symmetry property r(— t) = T(t) of the noise 
r(i) was utilized [see (|A1[) in Appendix |A"] . 



Eq. (|20j) . which may cause significant numerical error in 
the low-frequency regime, we can expressing the noise 
spectrum as a function of the AC component Q^f 2 (^) = 



T ft (uj) 



0.5 



(21) 



from Eq. (fig)) . 

The above formulas suggest that an identification 
scheme can be devised as shown in FigHJ For illustration, 
the formula (|19[) is used to identify the noise T(t) from 
the experimental data collected in NEC and RIKEN [sec 
Fig. EJb)] for testing coherent oscillations at the optimal 
point. The average charge number was measured with a 
time delay sweeping over a time period of T — 2900 ps 
with sampling time 8t = 9 ps. From the Laplace trans- 
form (at s = iui) of the oscillation curve, the character- 
istic frequency is read as A w 6.0 GHz from the peak 
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FIG. 4: Schematic diagram for noise identification in a super- 
conducting charge qubit. The average charge number Q(t) 
is measured and transformed into the frequency domain, and 
then used to calculate the noise spectrum via identification 
formulas (|20|) or (|2I[) at the optimal point. 



in Fig. [5ja). Then the spectrum of T(t) was calculated 
Eq. (|2l"j) as shown in Fig. 02b). 

This example illustrates how the identification strategy 
works with a simple experiment design without using ad- 
ditional control pulses or devices. The identification only 
uses the measurement data at the optimal point, but the 
samples should cover a sufficiently broad time interval. 
The identification error could come from many possible 
factors, including the measurement noise (not the noise 
coupled to the qubit) from the SET device. Moreover, 
the data collected from the experiment has a finite num- 
ber of data points and finite sampling time period, which 
restricts the precision and frequency regime of the iden- 
tification results. To improve the identification over a 
wider range, the data should be more densely sampled 
over a sufficiently long time period until the oscillation 
decays to approximately zero. 

Moreover, the linewidth of the coherent oscillation 
spectrum shown in Fig. @|a) is comparable to the width 
of Tft(w) in Fig. E^b), which seems to imply that the 
noise correlation time is comparable to the decay time. 
This could be due to the strong noise coupling (where 
the Born approximation, as well as the above derived 
identification formulas, are not valid), but it could also 
be from measurement errors (including the measurement 
noise and the finite number of data points). These factors 
should be further analyzed. 



IV. CONCLUSION 

We have presented a noise-identification approach that 
can improve the precision based on a frequency-domain 
model for the system's non-Markovian dynamics with less 
amount of measurement data. With applications to qubit 
systems, the Fermi's golden rule can be naturally derived 




FIG. 5: 



(color online) (a) The absolute value of the Laplace 

AC , 



transform of Q^%^p) ~ Q-n/iif) — 0.5, where Q n /2{t) is the av- 
erage charge number at the optimal point, (b) The identified 
noise spectrum Fft(i^) at the optimal point. 



under proper approximations from this model. The ad- 
vantage of this model is that a simple identification strat- 
egy can be obtained without aid of additional control 
pulses or devices. The obtained identification formula at 
the optimal point is illustrated with applications to su- 
perconducting qubits, which requires only measurements 
at the optimal point. 



It was also seen that the quality of the identification is 
limited due to the finite data points and possibly severe 
measurement noises in the experiment. Moreover, the ef- 
ficient identification of the noise term f2(t), which is non- 
ncgligiblc under ultra-low temperatures, is to be studied. 
In principle, the full expression (|C2[) derived in this paper 
can be used to extract the noise spectrum of f2(i) with 
data obtained under biased configurations. However, so- 
phisticated numerical algorithms are to be designed due 
to the complexity of the formulas. It is also possible to 
design other measurement schemes with which the iden- 
tification formulation can be simplified. These problems 
are to be explored in the future. 
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Appendix A: Properties of the noise correlation 
function 

According to the definition ([2]) of the noise correlation 
function, it is easy to examine the following symmetry 
properties: 

r(-r) = r(r), n(-r) = -n(r), $(-t) = #*(t). 

(Al) 

Owing to these symmetry properties, the Fourier trans- 
form of the symmetric function T(t) is a purely real- 
number valued function, while the Fourier transform of 
the antisymmetric function fl(r) is purely imaginary. 

In this paper, both the Laplace and Fourier transforms 
are involved to describe the noise spectrum. The major 
difference between them is that the Fourier transform is 
a two-sided integral in the time domain from r = — oo 
to r = oo, while the Laplace transform is one-sided (i.e., 
from r = 0~ to r = oo). Because $(r) is evaluated for 
both negative and positive times, the noise spectrum is 
generally defined as the Fourier transform of the corre- 
lation function $(r) pit [2o|. However, in the frequency- 
domain model (|7J) derived from the generalized master 
equation (JTJ), the Laplace transform has to be adopted 
because only the positive time-correlation function is in- 
volved. Another reason for using the Laplace transform 
is that the system is always prepared at some initial 
state in any realistic measurement, regarding to which 
the Fourier transform is physically not applicable. Nev- 
ertheless, owing to the symmetry properties, it is suffi- 
cient to use the positive-time branch, from which we can 
recover full noise correlation function. 

Therefore, the Laplace transform will be adopted cor- 
responding to the positive time-correlation function. To 
avoid confusion, the Fourier transform of a time- variant 
function, say will be denoted by $ft(w) with a 

subscript. $ft(w) is generally different from the Laplace 
transform evaluated at s — iu, but the symme- 

try property of the correlation functions guarantees the 
following relations: 

r FT H = 2ReT(iw), FT (o;) = 2Imft(iuj), (A2) 

which will be used to extract the temporal correlation 
function from the noise spectrum obtained via the inverse 
Fourier transform. 



Appendix B: Frequency domain Derivations in Sec. 
DTK 



Firstly, suppose that the eigenvalues of Cq are 
Xi, ■ ■ ■ ,x n , and that Cq is diagonalizcd by P, then the 
matrix s-function r(sl — Co) can be calculated as follows: 

r(sI-£ ) = p- 1 dia.g{T(s-x 1 ),--- ,T{s-x n )}P, 

and in a similar way ft(sl — Co) can be evaluated. 

With respect to the Hamiltonian (O, the system is 
initially prepared at the excited state corresponding to 
the Bloch vector v° = [1,0,0, —1] T . Correspondingly, 
the matrix representation Cq of Aer z is diagonalized as 

Co = P _1 diag{0, 0, iA, -iA}P, 



where 



P = 



10 
110 
-i i 
1 



(Bl) 



and hence 

r( S 7-£ ) = P" 1 diag{r( s ),r( s ),r(s-iA),r(s+iA)}p. 

When the qubit is released from the excited state, from 
Eq. ([5]) , the corresponding transition rate under the bias 
angle 6 can be expressed as 

a„ N (s) cot 2 8 + Ni(s) , , 

7i (s) = s n 2 6 ■ —^-4 9 ) { , (B2) 

m ; A,(s)cot 2 6> + L>i(s) 

where, by dropping the argument "s" for simplicity, we 
have 

N (s) = (s 2 + A 2 )r + + s(T 2 + +T 2 _) 

(ft - fl + )(AT + - sT-) - fl-(s 2 + sT+ + A 2 ) 

Ni(s) = (s 2 + s r + A 2 )(r+-ft_) 
d (s) = s [( s + r+) 2 + (r_ + A) 2 ] 
Dx( s ) = (s + r + )(s 2 + sr + a 2 ). 

Here, we use the following notations: 

r( s + iA) + r(s~iA) 
r+(«) = g 

= J8f [f(t)cosAt], (B3) 
T(s + iA) - r(s - iA) 
2i 

= J2f [7(t)sinAi], (B4) 

and in the same way are 0±(s) defined. These represent 
signals modulated by sinusoidal waves with the carrier 
frequency A. Such noise spectral functions shifted by the 
oscillation frequency ±A exhibit the interplay between 
the qubit and its environments. 

When the qubit is prepared in the ground state, the 
expression of 7f(s) is similar to Eq. (|B2[) except that the 
signs of all terms of fl(s) and f2±(s) are flipped. This 
shows that the presence of ft(t) in the noise correlation 
function causes the difference between the static average 
charge number and the corresponding transition rates. 



r_( s ) 
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Appendix C: Frequency domain Derivations in Sec. 



The extended Bloch vector corresponding to the ini- 
tial zero-charge state with Bloch vector being v° = 
[1, 0, 0, 1] T . We first diagonalize £q = Aag as 

C = P^diagjO, 0, iA, -iA}P e , 



where 



Pa 



10 

tan 8 — cot 8 — cot 6 

i esc 8 — icsc( 
11 1 



(CI) 



Consequently, we have 

r(s/-£ ) = P^Miag^^r^r^A^rXs+zA)}^, 



and from Eq. ([6]), the average charge Qe(t) for the bias 
angle 8 can be expressed as 



A JV ( s ) COS g + iVi(8) , . 



2s D (s) cot 2 + £>!(«)' 



where 



N (s) 
Ni(s) 
D (s) 
Di(s) 



(s + r+) [o + - n(*)] + (A + r_)n_ 
A( s + r + ), 

s [( s + r + ) 2 + (A + r_) 2 ], 
( s + r + )(s 2 + s r + A 2 ). 
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